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ABSTRACT 

We carry out ray tracing through five high resolution simulations of a galaxy cluster to 
study how its ability to produce giant gravitationally lensed arcs is influenced by the collision 
cross-section of its dark matter In three cases typical dark matter particles in the cluster core 
undergo between 1 and 100 collisions per Hubble time; two more explore the long ("colli- 
sionless") and short ("fluid") mean free path limits. We study the size and shape distributions 
of arcs and compute the cross-section for producing "extreme" arcs of various sizes. Even 
a few collisions per particle modify the core structure enough to destroy the cluster's ability 
to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to 
unrealistically large masses before it regains the ability to produce giant arcs. None of our 
models with self-interacting dark matter (except the "fluid" limit) is able to produce radial 
arcs; even the case with the smallest scattering cross-section must be scaled to the upper 
limit of observed cluster masses before it produces radial arcs. Apparently the elastic col- 
lision cross-section of dark matter in clusters must be very small, below 0.1 cm^g^', to be 
compatible with the observed ability of clusters to produce both radial arcs and giant arcs. 
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1 INTRODUCTION 

Several recent observations indicate a potential problem with the 
Cold Dark Matter (CDM) scenario of structure formation in that 
they seem to contradict the numerically simulated density struc- 
ture of CDM haloes. A first difficulty is posed by published rota- 
tion curves of dwarf galaxies, which often rise linearly from the 
centre to radii greater than '-^ 1 kpc. Such a rise implies a constant- 
density core f or the host dark-matter halo (Moore 1994 ; Flores & 



Primack 1994), and is in conflict with the central density cusps 
of the haloes found in W-body simulations, which are usually fit- 
ted by a double power-law profile with central logarithmic slope 
~ 1 ^ 1.5 (Navarro, Frenk & White 1997, hereafter NFW; Moore 
et al. 1999b). A second problem is related to the large number of 
low-mass subclumps found orbiting within simulated galaxy-size 
haloes; with a naive scaling to the observed properties of dwarf 
galaxies within the Local Group, the ab undance of these sub clumps 
appears to o high by a factor of 10 to 50 (|Moore et al. 1999at Klypin 
et al. 1999). 

A possible way to modify the CDM scenario to be compatible 
with these observ ations has recently been suggested by Spergel & 
Steinhardt (200C). They propose a model where the dark matter 
is self-interacting but dissipationless. They argue that if the dark 



matter particles have a sufficiently large cross-section for elastic 
scattering, the strength of the central density concentration and the 
abundance of orbiting subclumps can both be reduced. Collisions 
between dark matter particles would destroy the central cusps and 
lead to a gradual evaporation of small subclumps. 

Self-interacting dark matter would also affect the structure 
of galaxy clusters, making them more nearly spherical and re- 
ducing the number of subclumps. These effects have been con- 
firmed by several recent studies in which the halo properties in 
cosmological simulations of coUisionless cold dark matter and self- 
intera cting da rk matter have be en co mpared. For instance, Yoshida 
et al. ( |2000b| ) and Dave et al. (EOOll) found that intermediate scat- 
tering cross sections for dark matter particles produce haloes which 
are less centrally concentrated than in a coUisionless model and 
have smoother and more s pherica l cores. Similar results have also 
been obtained by Bi irkert ( |Z000| ), Firmani et al. (2000,2001) and 
Kochanek & White (200C ), although they find different time scales 



for the onset of core collapse, which transforms intermediate flat 
cores into cores with density profiles oc r^^. 

It is natural to ask how self-interacting dark matter would af- 
fect the strong gravitational lensing properties of galaxy clusters. 
Strong lensing is a very powerful tool to probe the distribution of 
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dark matter in cluster cores. From the location of multiple images 
of background sources lensed by a cluster one can obtain a de- 
tailed map of the ma ss distribution in the lens ( Tyson, Kochanski 
& deir Antonio 1998). The morphology of long arcs observed in 
t he cluster fields yields important information on the core density 
(Bartelmann et al. 1998). Finally, the location of radial arcs can put 



strong constraints on the size and compactness of the lens cores 
(e.g. Narayan & Bartelmann 1997). 

The goal of this paper is to evaluate whether self-interacting 
dark-matter models are compatible with the observed ability of 
galaxy clusters to produce strong lensing events. So me rece nt work 
has already addressed this issue. Miralda-Escude (2000) argued 



that the shallow cores produced by self-interacting dark matter 
could not agree with curre nt ob servations of lensing by clusters. 
Wyithe, Turner & Spergel ( |2000| ) studied simplified analytic mod- 
els of self-interacting dark matter haloes and found large collision 
cross-sections to be incompatible with multiple imaging by clus- 
ters. In this paper, we address the issue by using ray-tracing through 
high resolution simulations to compute lensing cross sections for 
arcs of various types. We study a set of simulations designed to 
isolate the effect of collisions between dark matter particles on the 
structure of an individual cluster. 

The plan of our paper is as follows. In Section ^ we describe 
the cluster models used in the simulations and present the numeri- 
cal method adopted to compute the lensing properties. We present 
and discuss the results of our analysis in Section ^. Finally, our 
conclusions are drawn in Section ^. 
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Figure 1. Surface density profiles of the five halo models at redshift z = 
0.278. The horizontal dashed line indicates the critical surface density for a 
lens at the same redshift and sources at redshift z = 2. Note that large arcs 
are formed at surface densities below the critical value because of the shear. 



2 NUMERICAL SIMULATIONS 



2.1 Cluster models 



The cluster models used in this study are those described in Yoshida 
et al. (2000a,b). They consist of a set of five different resimulations 
of the same cluster-sized halo, which is the second most massive 



object in the GIF-ACDM simulation of Kauffmann et al. ( |1999| ). 
All simulations were performed with the GADGET code (Springel, 
Yoshida & White 2000), in a flat model universe, with a matter 
density parameter i2,n = 0.3 and a cosmological constant of flyv = 
0.7. The Hubble constant is assumed to be Hq = 70km s^' Mpc^' . 
The CDM power spectrum was normalised so that the rms matter 
density fluctuations in spheres of radius r = 8/z^' Mpc is Og = 0.9 
(h is the value of Hq in units of 100 km s^' Mpc^'). 

The first cluster of the series (hereafter called the SI model) 
is just a resimulation at higher resolution of the original GIF- 
ACDM cluster, and was produced assuming standard collisionless 
dark matter. Three other simulations (hereafter called W-models) 
were carried out introducing elastic scattering between CDM parti- 
cles. This w as acc omplished by the Monte-Carlo method proposed 
by Burkert (20001, slightly modified by considering the pairwise 



velocity difference in evaluating the scattering probability instead 
of the one-point velocity dispersion. The three simulations dif- 
fer in the value of the scattering cross section per unit mass o*: 
O.lcm^g^^, l.Ocm^g^' and lO.Ocm^g^' were chosen. Hereafter, 
we will refer to these models as SlWa, SlWb an d SlW c, respec- 
tively, following the notation of Yoshida et al. (2000b). Finally, 



Yoshida et al. (2000a) also considered the extreme case in which 



dark matter particles are as strongly interacting as in the "fluid" 
limit (see also Moore et al. 2000). In this last simulation, the time 
evolution, starting from the same initial conditions as in the colli- 
sionless original simulation, was followed solving the fluid equa- 



tions with the SPH technique. In the following, this model will be 
referred to as SIF. 

The resolution achieved in these simulations is quite 
high. Most of them employ 0.5 x 10^ particles in the cen- 
tral, high-resolution region, where the particle mass is nip = 
0.68 X 1O'"/?^'M0. The gravitational softening length is set to 
20/i^'kpc. All haloes have a similar final virial mass of 7.4 x 
lOl4/;-iM0. 

For our lensing analysis, we picked the simulation snapshots at 
redshift z = 0.278 because this is close to the most efficient redshift 
for strong gravitational lensing (e.g. Bartelmann et al. 1998). 

The two-dimensional density profiles of the five cluster mod- 
els at this redshift are shown in Fig. |^. The general trends in the 
profiles follow our expectations. As the scattering cross section of 
the dark matter particles increases, the cluster density profile in the 
inner region of the halo flattens. Note that the central density for 
the model S IWb is smaller by more than one order of magnitude 
than that of the model S 1 . 

The halo shape also depends on a*. The principal axis ratios of 
the five models, determined from the inertial tensors of the matter 
at densities exceeding 100 times the critical value, are listed in the 
third column of Table |l[ The values clearly show that the haloes 
are more spherically symmetric when the frequency of collisions 
between the dark matter particles increases. In the fourth column 
of Table [l[ we also report the core radius rc of the three haloes with 
self-interacting dark matter, which we define as the distance from 
the cluster centre where the surface density drops below half its 
central value. For models SlWb and SlWc we typically find core 
radii larger than 200/;^' kpc. 
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Table 1. The main properties of the cluster haloes at redshift z = 0.278. 
Column 1: halo model; column 2: scattering cross section per unit mass c,; 
column 3: axial ratios; column 4: core radius r,., defined as the clustercentric 
distance where the surface density falls below half its central value. 



Radial and tangential critical lines are located where the con- 
ditions 
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SI 




1:0.58:0.54 
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SlWa 
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1:0.63:0.60 


110 


SlWb 


1.0 


1:0.67:0.65 


245 


SlWc 


10.0 


1:0.86:0.85 


400 


SIF 




1:0.95:0.81 


35 



2.2 Lensing properties of the clusters 

Our method to study the lensing properties o f gala xy clusters has 
been described in detail in Meneghetti et al. (2000). Compared to 
that paper, only few parameters are changed here; therefore, we 
only give a brief description here and refer the reader to the original 
paper for further details. For a general introduction to the th eory o f 
gravitational lensing, se e, e.g . Schneider, Ehlers & Falco (1992), 



Narayan & Bartelmann (1997) and references therein. 

Starting from a cluster simulation, we construct cluster lenses 
as follows. We centre the cluster in a cube of 3/i^' Mpc side length 
and obtain the three-dimensional density field p by interpolating 
the mass density within a regular grid of 256-' cells, using the 
Triangular Shape Cloud method (TSC; see Hockney & Eastwood 
1988). We then produce three different surface-density fields Z per 
cluster by projecting p along each of the three coordinate axes. 
This yields three lens planes per cluster, which we consider in- 
dependent cluster models for the present purpose. As long as the 
TSC smoothing kernel is sampled by sufficiently many particles, 
smoothing and projecting can be interchanged, so that smoothing 
the three-dimensional density prior to projecting it is equivalent to 
projecting first and then smoothing the two-dimensional density. 
We checked and confirmed that changing the grid resolution up by 
a factor of two yields surface-mass densities with indistinguishable 
lensing properties. 

The surface density fields are then scaled by the critical sur- 
face mass density for lensing, which depends on the cosmological 
parameters and on the lens and source redshifts. We recall that the 
redshift for the lensing cluster is zl = 0.278, and we put all sources 
to zs = 2. This choice is motivated by the fact that galaxy clusters 
at redshifts between 0.2 and 0.4 are most efficient strong gravita- 
tional lenses for sources at redshifts z 1. Although real sources 
are distributed in redshift, putting all of them at a single redshift is 
admissible because the critical surface density changes very little 
with source redshift if the lens redshift is substantially smaller, as 
is the case here. 

Scaling with the critical surface density yields three two- 
dimensional convergence fields K for each cluster. We then propa- 
gate a bundle of 2048 x 2048 light rays through the central quarter 
of each of these fields. Their deflection angles, directly obtained 
from the convergence, are used to compute the shear y and the ele- 
ments of the Jacobian matrix A, which describes the local properties 
of the lens mapping. It is symmetric and can thus be diagonalised. 
Its two eigenvalues are 



= and A.1- = 



(2) 



are satisfied, respectively. The corresponding caustics in the source 
plane, close to which sources are imaged as large arcs, are obtained 
applying the lens equation to the critical curves. All lens properties 
are computed on grids with an angular resolution of 0.19" on the 
lens plane, so that lensed images are properly resolved. 

The sources are initially distributed on a regular grid in the 
source plane. Their spatial density is iteratively increased near 
caustic curves. Placing a larger number of sources where the lens 
strength is highest increases the probability of producing long arcs 
and thus the numerical efficiency of the method. In the following 
statistical analysis, it is then necessary to compensate for this arti- 
ficial increase by assigning to each image a statistical weight pro- 
portional to the inverse of the resolution of the grid on which its 
source was placed. 

Sources are assumed to be elliptical, with axis ratios randomly 
drawn from the interval [0.5, 1], and area equal to that of a circle of 
1" diameter. We checked that changing the average source size does 
not affect our final results. 

For the classification of imag es we follow the technique in- 
troduced by Barte lmann & Weiss (1994), which was also used in 
Meneghetti et al. ([2000[). 



Xt = 1 



- Y and = 1 — K + y . 



(1) 



3 RESULTS 

3.1 Distributions of image properties 

In this section, we present and compare the statistical properties 
of the images of background sources lensed by the models of our 
cluster sample. For each image, we compute the length / (defined 
as the maximum length of the circular segment passing through 
the image), the width w (found by fitting the image with several 
geometrical figures; for further details, see Meneghetti et al. 2000), 
and the length-to- width ratio l/w. For the following analysis, we 
combine the results of two sets of lensing simulations, in which we 
used the same lenses but changed the distribution of the background 
sources by drawing different orientations and axial ratios. With this 
procedure we confirm the robustness of the results and improve 
the statistics. Searching for large arcs, we can neglect the images of 
unlensed or even weakly lensed sources. Therefore, we can exclude 
those images from the analysis whose area is smaller than two times 
the area of their unlensed sources, which also avoids uncertainties 
in the automatic arc classification and parameter estimation. The 
number of images in the final analysed samples ranges between 
1 1,252 for model S IWc and 24,695 for model SIR 

We plot the distributions of arc properties in Fig. ^ The his- 
tograms show the number of arcs observed in the field of the five 
cluster models as a function of their length (top panels), width 
(middle panels) and length-to-width ratio (bottom panels). Each 
column refers to a different model. From left to right; the pure col- 
lisionless case SI, the three models SlWa, SlWb and SlWc (with 
cross section of 0.1, 1.0, and lO.Ocm^g^', respectively), and the 
"fluid" limit model SIF. To facilitate the comparison, we show in 
Fig. I the 0, 1, 10, 50, 90, 99, and 100 percentiles of the distribu- 
tions plotted in Fig. ^ 

Regarding arc lengths, we find that model SI produces a 
higher number of long arcs than the collisional W- models. More- 
over, as shown in the top left panel of Fig. ^ the highest percentiles 
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Figure 2. Histograms of the image lengths (top panels), widths (middle panels) and length-to-width ratios (bottom panels). Images, whose area is smaller than 
two times the area of the unlensed source, are not considered here. Each column refers to a different model; from left to right: the pure collisionless case SI, 
the three models SlWa, SlWb and SlWc, and the "fluid" hmit model SIR 



(i.e. the largest arc lengths in the sample), change significantly go- 
ing from Si and SlWa to SlWb and SlWc. This is easily under- 
stood in terms of the different density profiles of the models in 
the inner regions of the cluster. S 1 has a steeper profile, and so it 
reaches a higher central surface density than models SlWa, SlWb 
and SlWc. The tangential critical curve is located where = 0, 
i.e. where K = 1 — y, and K is just proportional to the lens surface 
density. Therefore, the curve shrinks as the surface density pro- 
file flattens. Consequently, the probability of producing long arcs 
is highest for model SI, which has the longest tangential critical 
curve. 

Models SlWb and SlWc do not produce arcs longer that ~ 3", 
hence these models are unable to strongly distort background im- 
ages. Their central surface density is 55% and ^ 50% of the 
critical value, respectively, and the shear is reduced by the reduc- 
tion of subhaloes and the increased axial symmetry caused by dark- 
matter interaction. Consequently, the combination of convergence 
and shear never succeeds to produce critical lines, regardless of the 
projection considered. At the opposite end, the fluid limit case SIF 
yields a number of long arcs larger than SI because the extremely 
high density in the centre of this cluster pushes the tangential criti- 
cal curve outward and extends it. 

The distribution of image widths is also sensitive to the size of 
the interaction cross section. The zero-percentile (giving the width 



of the thinnest arcs in the sample) moves to higher widths as Ov, 
increases, indicating that SI produce s thin ner arcs than the W - 
models. As demonstrated by Kovner (1989) and Hammer (1991), 



the radial extension of a tangential arc is magnified by a factor of 
order /Jr ~ [2(1 — k)]^' relative to its original size. Therefore, if 
K ^ 0.5 at the tangential critical line, tangential arcs are demag- 
nified in width. This condition can be achieved only if the central 
density of the lens is particularly high with a steep surface-density 
profile, as is the case for model SIF, or if the lens has an aspher- 
ical mass di^stributioi^^o^^^ronoun^ & 
Weiss 1994T TBartehnanner^L'T995^ !'7i^^ lens features 

both contribute to pushing the critical curves outward, where K is 
lower. As the steepness of the density profile and the predominance 
of substructures both decrease going from model SI to the W- 
models, the observed increase of the zero-percentile shown in the 
second panel of Fig. ^is consistent with our expectations. The lack 
of substructures in model SIF is compensated by the very dense 
core, providing this cluster with the capability of producing very 
thin arcs. 

From the middle panels in Fig. ^ and the middle panel of Fig. ^ 
(100-percentile), it can also be seen that the number of thick arcs 
decreases as the interaction cross section grows. The width dis- 
tribution for SlWb and SlWc reflects the fact that these two mod- 
els produce only weak distortions on the images of the background 
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Figure 3. Percentiles of the arc length, width, and length-to-width ratio 
distributions are shown at the levels of 0% (open triangles and dashed lines), 
1% (filled triangles and dotted lines), 10% (filled pentagons and dotted- 
dashed lines), 50% (filled squares and solid lines), 90% (open circles and 
long-dashed lines), 99% (filled circles and dotted lines) and 100% (heavy 
open circles and dashed lines). The figure shows the lengths (or widths or 
length-to-width ratios) below which fall 0, 1, 10, 50, 90, 99, and 100 per 
cent of the simulated arcs for the different cluster models. For instance, it 
can be read off the top panel that all arcs produced by model S 1 are shorter 
than 25". 



sources, that is, both image lengths and widths are only shghtly 
changed with respect to their intrinsic values. On the other hand, 
the higher central surface density of the SI, SlWa and SIF models 
makes them critical for strong lensing. 

Finally, and consistently with the previous discussion on arc 
lengths and widths, the distributions of the image length-to-width 
ratios show that, when the dark matter particles interact more fre- 
quently and the density profile becomes shallower, the number of 
arcs with large length-to-width ratio decreases. That is not the case 



for model SIF, where the density profile is even steeper than in 
model SI, causing the creation of many very long and thin arcs. 



3.2 Lensing cross sections for long and tliin arcs 

We now compute absolute lensing cross section of the clusters for 
long and thin arcs. Arc cross sections for a property Q are defined 
as the area in the source plane within which a source has to lie in 
order to be imaged as an arc with property Q. We consider here the 
cross sections for arcs with length-to width ratio larger than three 
different thresholds, namely: 



Qi 

Qi 

23 



(//w) > 5 
(//w) > 7 
(//w) > 10 



In particular, the i mages which satisfy co ndition are commonly 
called giant arcs (Wu Sl Hammer 1993). We are particularly in- 
terested in arcs with large length-to-width ratio because they have 
been observed in many galaxy clusters, and their statistics has been 
used in previous studies in an attempt to constrain cosmological 
parameters (Bartelmann et al. 1998 and references therein). 

Each cluster model provides six estimates of its strong-lensing 
cross sections, two for each projection since we performed simula- 
tions with two source distributions. We average all measurements 
in order to obtain the mean cross sections for each cluster model. 
The lensing cross sections of the five models are listed in Tab. |^. 
Models SlWb and SlWc have vanishing cross sections even for 
arcs satisfying condition Q\ , given that, as shown in Fig. ^, they do 
not produce arcs with length-to-width ratio larger than ~ 3.5. 

We find that model SI has cross sections larger than model 
SlWa. In particular, the cross sections Oq, and ctq, of the first 
model are more than one order of magnitude larger than those of 
the second. This feature suggests that a galaxy cluster consisting 
even of very weakly interacting dark matter should produce at least 
ten times fewer giant arcs than one made of collisionless particles. 
Compared to the other models, SIF has very large cross sections, 
in agreement with the results shown earlier. 

As a consequence of the lower curvature in the effective lens- 
ing potential, arcs produced by flatter profiles tend to be more 
strongly magnified. However, the corresponding change in the 
magnification bias is expected to have a very small effect on the 
previous results for two reasons. First, such arcs would then be 
very thick, while the observed ones are generally thin and some- 
times unresolved in the radial direction even when observed from 
space. Second, the magnification bias depends on the colour selec- 
tion, because the number-count slope does, and unless only very 
blue or very red objects are selected, the bias is negligible for arcs. 

A comparison of the cross sections listed in Tab. ^ for the col- 
lisionless simulation (SI) shows that the cluster almost exactly re- 
produces the average cross section found in the cluster sample used 
by Bartelmann et al. (1998) for a lens redshift of 0.28. In addition, 
we checked whether different realisations of the same cluster had 
comparable large-arc cross sections. To do so, we took the clus- 
ter mass distributions at different output redshifts {z = 0.20, 0.13 
and 0.06) and computed their lensing properties when placed at a 
redshift of 0.28. For the given cosmological parameters, the linear 
growth factor changes by less than 20% between redshifts 0.28 and 
0.06, so that the effect of cluster growth is unimportant. We find 
only a modest growth in the large-arc cross sections when using the 
cluster mass distributions at redshifts 0.20 and 0.13. Between red- 
shift 0.13 and 0.06, however, the cluster undergoes a merger event 
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Table 2. Lensing cross sections of the five halo models. Column 1: halo model; column 2: cross section c^i for ^^s with l/w > 5; column 3: cross section 
0g, for arcs with l/w > 7; column 4: cross section Cq^ for arcs with l/w > 10; column 5: cross section CQf, for radial arcs. All cross sections, given in units 
of (Mpc//))^, are averaged over the three projections of each cluster and two source-plane realisations for each projection. 



Model Gq, Cp, Gqi^ 

SI 6.93x10-3 4.82x10-3 2.20 x lO"' 8.13x10-" 

SlWa 1.46x10--'' 2.37x10-" 1.26x10-" 

SlWb 

SlWc 

SIF 3.83x10-2 2.36x10-2 1.14 xUr^ 1.02 xlO-^ 



with a clump of 10 Mr., and subsequent relaxation, which in- 
creases its central density and therefore also the large-arc cross 
section. The effect of the merger on the lensing cross section is 
larger for the coUisional than for the coUisionless model. The in- 
falling clump increases surface mass density and shear. Since the 
surface density profile of the coUisional cluster is flatter than that 
of the coUisionless cluster, this increase in lensing efficiency af- 
fects a larger region in the cluster centre. While the ratio between 
the large-arc cross sections between the SI and SlWa models is 
roughly an order of magnitude during quiescent periods, it is re- 
duced to approximately five during the merger event. For similar 
reasons, it is important to mention that the presence of a dominant 
central galaxy could affect our different models in different ways, 
being more effective on top of the flatter models. As a consequence, 
the differences between coUisional and coUisionless cases could be 
reduced. A detailed study of this effect is in progress. 

Nonetheless, we can conclude (1) that the chosen cluster is 
typical for clusters in the ACDM cosmology, and (2) that the large- 
arc cross section for the mildly coUisional model SlWa falls sub- 
stantially below that of the coUisionless model S 1 for different re- 
alisations of the same cluster. 
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3.3 Radial arcs 

In this subsection, we analyse the capability of our models to pro- 
duce radial arcs. Such arcs have been observed in several galaxy 
clusters. For insta nce, two radial arcs exist in the core of MS 0440 
( [Oioia et al. 199^ ). From X-ray observations, the mass of this clus- 
ter within ~ 350/;-' kpc was estimated to be 0.65 x lO'"* /?-' Mq, 
which is similar to that of the cluster we analyse here. 

The appearance of radial arcs and their position within the 
lensing cluster depend on both the slope of the projected mass pro- 
file and on the central density of the lens. The surface density must 
surpass a critical value in order for the radial critical curve to appear 
where the eigenvalue ^ = 1— K + yof the Jacobian matrix is zero. 
Moreover, the steeper the density profile is, the closer to the cen- 
tre t he radial arcs tend to move (W illiams, Navarro & Bartelmann 
1999). The observed location of radial arcs close to the cluster mass 
centre s suggests that the core radii of the host clusters m ust be quite 
small dMellier, Fort & Kneib 1993t |Smail et al. 1996|; H ammer et 
al. 1997; pioia et al. 1998| ). Due to the different slopes of the den- 
sity profiles of the cluster models, we thus expect to find quite large 
differences in their capacity to form radial arcs (see e.g. Molikawa 
& Hattori 2000). 

The radial magnification fj^ at each point of the lens plane is 
given by the inverse of the radial eigenvalue of the Jacobian lensing 
matrix, /j^ = l/A^. From the convergence and the shear, computed 
as explained in Sect. 2.2, we obtain /jj- at each of the 2048^ grid 
points of the lens plane. We show in Fig. ]A how these radial magni- 



Figure 4. Histograms of the radial magnifications fj^ at each point on the 
lens plane for the five cluster models: SI (solid line), SlWa (dotted line), 
SlWb (dashed fine), SlWc (long-dashed line) and SIF (dot-dashed line). 



fications are distributed for the five cluster models. Each histogram 
combines contributions from all three projections of each cluster. 
As can be noted, values of exceeding 10-' are reached only by 
models SI and SIF, while models SlWa, SlWb and SlWc can 
only produce weak radial distortions. 

Using the tangential and radial eigenvalues of the Jacobian 
matrix, we can also identify the radial arcs from the complete sam- 
ple of distorted images, and this enables us to compute the lensing 
cross sections for radial arcs. For their automatic detection, we as- 
sumed that radial arcs must lie on the radial critical curve or very 
close to it. Therefore, we searched those images containing at least 
one pixel for which the radial eigenvalue is < 0.15 and the tangen- 
tial eigenvalue is > 0.6. These values were chosen in order to avoid 
misclassifications of radial arcs. 

The cross sections for radial arcs. Or, averaged over the six 
simulations performed for each cluster, are listed in the fifth col- 
umn of Table ^ We found that only models SI and SIF are able 
to produce radial arcs, while the W-models are not. The results for 
model SlWa are particularly interesting because this cluster seems 
to have a sufficiently high surface density to produce giant arcs, but 
not enough to radially distort the background galaxies. 
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Figure 5. Lensing cross sections of cluster model SlWb as a function of 
the rescaling factor /. The three panels show the cross sections for arcs 
with l/w > 5 {<SQ^, top frame), l/w>l (<5q^, middle frame) and l/w > 10 
(Ogj, bottom frame). The dashed and dotted lines in each panel indicate 
the lensing cross sections of cluster models SI and SlWa, respectively. The 
cross sections are given in units of (Mpc//i)^. 



3.4 Rescaling the clusters 

Our results show that models S 1 Wb and S 1 Wc cannot produce long 
and thin arcs at all. Although the mass of these clusters is already 
fairly large, we have checked whether an even more massive cluster 
with density profile similar to that of SlWb could efficiently form 
this kind of arcs. To this end, we rescaled the cluster mass by a 
variable factor / > 1 . 

The virial radius i?vir of a halo scales as i?vir °^ M^l^ , where 
M is the virial mass of the halo. To obtain a halo of mass f x M 
which is dynamically stable, we also need to rescale distances by a 
factor /'/^. This means that, while the three-dimensional density p 
remains fixed, the halo surface density is enhanced by a factor /'''"'. 
For this reason we expect that, increasing its mass, cluster SlWb 
will eventually become able to produce strong lensing effects. 

We therefore perform a new set of lensing simulations, choos- 
ing mass scaling factors / equal to 3, 5, 7, and 10. We plot in Fig. ^ 
the cross sections of these rescaled clusters, as a function of the 
rescaling factor /. In this plot, upper, middle and lower panels refer 
to the lensing cross section for length-to- width ratios larger than 5, 
7, and 10, respectively. The dashed and dotted horizontal lines cor- 
respond to the equivalent cross sections for models SI and SlWa. 

We can infer from this figure that a cluster with the same den- 
sity profile as our model SlWb must have a mass of at least ~ 5.5 
(~ 4.5) times larger than the mass of our present cluster, in order to 
achieve the same strong lensing efficiency of model SI (SlWa). 

In the same spirit, we also rescale the mass of cluster SlWa, 
to determine when it starts producing radial arcs. We multiply the 
cluster mass by a factor of two and found that only one out of three 
projections of the rescaled cluster was efficient in forming radial 
arcs. However, the cross section of the rescaled cluster, Gr = 2.68 x 
10"^ is still one order of magnitude lower than that of model S 1 . 



Finally let us comment that our method to rescale the clusters is 
approximate. In fact, rescaling the mass should in principle affect 
the collision rate. Yoshida et al. (2000b) estimated the number of 
collisions between dark matter particles per Hubble time as A'coii = 
pOicV, where p is the three-dimensional density and V is the particle 
velocity. Since V <>= M'/^, the number of collisions would increase 
by a factor /'/^ when the cluster mass is rescaled by a factor /. 
Since our rescaling technique keeps N^oll constant while the cluster 
mass is increased by a factor /, the result practically corresponds 
to an effective reduced by a factor /'/^. This means that we are 
even over-estimating the lensing efficiency of the /-scaled clusters 
and our previous conclusions can be considered conservative. 



4 SUMMARY AND CONCLUSIONS 

Self-interacting dark matter was suggested as a solution to the 
problems (1) that dwarf galaxies are observed to have flat cores, 
while they should have density cusps according to A'-body simu- 
lations, and (2) that simulations tend to produce more substructure 
in galaxies than there seems to be. In essence, the dark-matter self- 
interaction flattens density cusps, increases radial symmetry, and 
damps substructure; all of which are 'desired' effects. 

However, it was shown before that asymmetries, steep density 
profiles, and existence of substructures are essential for the ability 
of galaxy clusters to produce pronounced strong-lensing features 
like large arcs. This is due to several reasons. First, lenses need to 
exceed the critical surface mass density in order to form critical 
curves, which are mandatory for large arcs. Second, asymmetries 
in lenses increase the gravitational tidal field, or shear, which helps 
to make them critical for multiple imaging. Asymmetric clusters 
can be critical at lower surface mass densities than symmetric clus- 
ters. Third, observed arcs are generally thin, i.e. little magnified or 
even demagnified in the radial direction. For this to be the case, 
the surface mass density needs to be lower than about half the crit- 
ical surface density at the location of the arcs. Together with the 
requirement that the surface density be supercritical in lens centres, 
this implies that density profiles have to be steep. 

In other words, the introduction of self-interaction between 
dark-matter particles is expected to have a pronounced, potentially 
'undesirable' impact on the ability of clusters to produce large arcs: 
it is exactly the desired effect of the self-interaction, namely to 
make clusters flatter, rounder, and smoother, that on the other hand 
threatens to destroy their strong-lensing capabilities. In this paper, 
we have investigated this effect on the second-most massive cluster 
in a cosmological ACDM simulation, and produced the following 
results: 

• The introduction of self-interaction even with a small cross 
section of = 0. 1 cm^ g^ ' reduces the cluster cross section for the 
production of large arcs by about an order of magnitude compared 
to the collisionless case. 

• Upon further increase of the interaction cross section, the clus- 
ter becomes entirely uncritical and cannot produce arcs any more, 
despite the otherwise favourable conditions in terms of cluster mass 
and redshift. 

• The cross section for radial arcs vanishes even for the smallest 
interaction cross section considered. 

• The cluster with an interaction cross section of = 
l.Ocm-^g^^ would require about five times more mass to reach 
the cross section for large tangential arcs of the cluster with o^, = 
0. 1 cm^ g^^ , and about 6 — 8 times more mass to restore the cross 
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section of the collisionless cluster. Due to the steepness of the clus- 
ter mass function, the spatial number density of such clusters at the 
redshift considered is more than six orders of magnitude smaller 
than that of the original cluster 

• The strong-lensing cross section of the cluster model inves- 
tigated here is almost identical to the average cross section found 
in a sample of clusters simulated earlier in the ACDM cosmology 
(Bartelmann et al. 1998). Also, the cross sections obtained do not 
change significantly when the cluster model is taken at different 
output times and placed at the original lens redshift of z = 0.28. 

At a length-to- width ratio of //w > 10, the abundance of ob- 
served arcs brighter than R = 22.5 is ~ 0.2 — 0.3 per cluster with 
X-ray luminosity > 2 x 10**ergs^' . Extrapolated to the whole sky, 
there should be ~ 1500 — 2300 such arcs. Numerically simulated 
clusters in a (cluster-normalised) ACDM model fall short by about 
an order of magnitude of producing this number of arcs (Bartel- 
mann et al. 1998). It may be possible that massive central cluster 
galaxies can increase the arc optical depth of ACDM clusters suf- 
ficiently to reconcile it with observations (Williams et al. 1999). 
However, a further reduction of the strong-lensing cross section 
by an order of magnitude due to even mild dark-matter interaction 
seems problematic in view of the observed arc abundance. 

Although the inability to produce radial arcs of the clusters 
with interacting dark matter is a potentially important piece of in- 
formation, it is currently impossible to draw any firm conclusions 
from comparisons with data. Radial arcs have so far been reported 
in five galaxy clusters (MS 2137, Fort et al. 1992; A 370, Smail et 
al. 1996; MS 0440, Gioia et al. 1998; AC 1 14, Natarajan et al. 1998; 
A 383, Smith et al. 2000), all of which are located near bright galax- 
ies in the cluster centre. It is currently unclear what exactly the 
influence of the galaxies on the occurrence of the radial arc is, and 
how the presence of luminous galaxies near cluster centres prevents 
the detection of more radial arcs. A statistical comparison of pre- 
dicted and observed radial-arc numbers therefore appears prema- 
ture. In principle, the ratio of the clustercentric distances of radial 
and tangential arcs could be used to constrain the cluster density 
profile, but in practice at least redshift estimates of the arc sources 
would have to be available. 

Of course, our simulations are insufficient for statistically 
sound statements on the total arc cross section of a given cluster 
population. However, given the massive impact of dark-matter in- 
teraction even with the smallest cross section considered, we con- 
clude that the dark-matter self-interaction hypothesis may be in se- 
vere conflict with the abundance of large arcs unless the interaction 
cross section is very small. 
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